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Abstract 

A new oscillator model with different form of the non- 
minimal substitution within the framework of the Duffin- 
Kemmer-Petiau equation is offered. The model possesses 
exact solutions and a discrete spectrum of high degeneracy. 
The distinctive property of the proposed model is the lack of 
the spin-orbit interaction, being typical for other relativistic 
models with the non-minimal substitution, and the different 
value of the zero-point energy in comparison with that for the 
Duffin-Kemmer-Petiau oscillator described in the literature. 
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1 Introduction 



The harmonic oscihator is one of the most basic and useful solvable examples in the 
non-relativistic quantum mechanics. Recently, the harmonic-oscillator potential for 
relativistic spin-1/2 particles has received considerable attention by many groups. 
Besides the academic importance of the exact solutions, the relativistic oscillator 
may be an analytical basis to deal with more realistic interactions in several areas, 
namely, in nuclear and particle physics. In addition to the exact solutions of the 
Dirac equation with confining scalar- vector potentials El 0] , since the paper 
by Moshinsky and Szczepaniak the interest has been extended to the oscillator 
model constructed by means of a non-minimal coupling scheme. 

This model, introduced formerly in the works [7| HI [5] , has the Hamiltonian 
which is linear in both momenta and coordinates, conserves angular momentum, is 
exactly solvable and its eigenspectrum is highly degenerate. The system obtained, 
being a harmonic oscillator with a spin-orbit coupling in the non-relativistic limit, 
is referred to as the Dirac oscillator (for further references see Ref. 10). 

Analogous models have been constructed within the framework of the Klein- 
Gordon equation [111 112) . and the spin-0 and spin-1 Duffin-Kemmer-Petiau (DKP) 
formalism, too. In the latter case, there are two different approaches. The first 
approach jl3[ 114) . the PSUSY oscillator, employs ideas of the parasupersymmetric 
quantum mechanics ^S] whereas the second one the DKP oscillator, exploits 
the technique developed for the Dirac oscillator [5]. However, the PSUSY oscillator 
has been proved to be equivalent to the DKP oscillator supplemented with the 
additional constraint [T7] . 

The eigenstates of the spin-1 DKP oscillator subdivide into the states with the 
natural and unnatural parities, whereas for the PSUSY oscillator only the natural- 
parity states become permissible due to the above-mentioned constraint. It should 
be pointed out that the energy eigenvalues of unnatural-parity states of the spin- 

1 DKP oscillator show the behavior which is unusual for oscillator models. These 
energy eigenvalues are not equidistant in contrast to those of the spin-0 and spin-1/2 
oscillators and the spectrum of the spin-1 oscillator with the natural-parity states. 

The goal of the present work is to study a more general type of the non-minimal 
substitution allowed within the framework of the DKP equation and to show that 
there is an alternative model of the spin-1 DKP oscillator whose energy spectrum 
is equidistant for the states of both the natural and unnatural parity. 

2 Generalization of the Model of the DufRn-Kemmer- 
Petiau Oscillator 

In this section we sketch out the technique of the non-minimal substitution and 
consider the generalization of the model of the DKP oscillator. 

Let us recall that the Duffin-Kemmer-Petiau equation |18[ I19L OH] describing a 
free scalar or vector boson of the nonzero mass m can be written as 

ih(3°^ = {cP-p + mc')^, (1) 

where the matrices /3'^ (/i = 0, 1, 2, 3) satisfy the Kemmer trilinear algebra 

P'^P'^P^ + P^P^p^' = gf'^P^ + g^'^P" (2) 

and the convention for the metric tensor is g'^'^ = diag{l, —1, —1, —1). 

Actually, there exist different representations for this algebra ^T]. However, the 
most used among them are the following three very well known cases: the trivial 
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one with (3^^ = 0, a, five-dimensional representation describing the spin-0 particles, 
and a ten-dimensional one for the spin-1 bosons. In the last case, the matrices /?'^ 
in the notation of the paper |16j are given by 
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where Si are the usual 3x3 spin-1 matrices and 
= (0 0) ei = (l 0) 62 = (0 



1 ) 63 - ( 1 ) (4) 



while and I stand for the 3x3 zero and identity matrices, respectively. 
Inserting the non-minimal substitution 

p ^ p — imLurjQr 

into we obtain the equation for the DKP oscillator j^H] 



ih(3' 



0^ 

dt 



[c/3 • {p — imujr]Qr) + mc^]^ 



(5) 



(6) 



where tjq = 2(/3'^)^ — 1 and u is the oscillator frequency. 

We attempt to generalize this model through a modification of the non-minimal 
substitution. For this purpose we extend the coupling term to the more general 
form 

p p ~ imojMr, (7) 

where M is some 10 x 10 matrix consisted of combinations of matrices j3^. 

Obviously, not every matrix M will lead us to an oscillator model. As it was 
noticed in Ref. 12, in models of both the spin-0 and spin- 1/2 oscillators, a matrix 
involving a potential in the non-minimal substitution anticommutes with the ma- 
trix structure of the momentum part of the wave equation. We suppose that this 
property remains valid as well for a model of the spin-1 oscillator. In our case, the 
such matrix is M, so that the following condition must be held 



f3M + M(3 = 0. 



(8) 



For the Kemmer algebra (0), there exist four matrices anticommuting with /3. 
Hence, M may be chosen as their linear combination 

M = ai?/o + a2% + as?/ + a4'7o?75?7i (9) 
where are numerical coefficients and the matrices 775 and ry are given by 
= -(2(/3Y-l)(2(/3^)' + l)(2(/32)2 + l)(2(/33)2 + i), 
77 = 2/3^/3^-5. (10) 

Exploiting the spin-1 representation for we have 

V = -V5, VoV5V^-Vo, (11) 
and the final expression for the matrix M reduces to 

M ^birjo + b2V5- (12) 
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Thus, the DKP equation with the non-minimal substitution Q which includes 
the matrix M in the form (|12|) is written as 

dtp 

ihP°— = [cP ■ (p- imLu{birio + hvs)^) + ■mc^]ip. (13) 
at 

Next, the block structure (|3J) of the matrices implies that the ten-component 
wave function of a stationary state may be decomposed into 

V'(r,t) = ( i(p(r) F(f) G(f) K{r) f exp {-iEt/H) (14) 

and the equation 11311 gives rise to the following set of time-independent equations 

mc^(p — ic{p — i{bi + h2)mwr) ■ G, 

mc^F = EG - c{p + i{bi - b2)mLur} X K, 

mc^G = EF + ic{p + i{bi + b2)mujf)(p, 

m(? K = —c{p-^ i(bi — b2)mijjr) x F. (15) 

Upon eliminating 1^(^)7 G{r) and K{r) this set becomes equivalent to the only one 
equation for the 3-component spinor F{r) 

{E^ - m2c4)F = [c^f + (61 - b2)^m^u;^f^) 

-(36i - b2)hujmc^ - 2biUjmc^L ■ s\F + 46i52C^m^w^f(f • F) 

261 w - ^ 

—{p + i{bi + b2)mLur)[{p- i{bi - b2)mu;f) ■ {2h + L ■ sjF], (16) 

m 

where L is the orbital angular momentum and s is the 3x3 spin-1 operator. 

The derived equation, apart from the last term, is that for a harmonic oscillator. 
However, as c ^ 00 this last term becomes negligible and in the non-relativistic limit 
the described system corresponds to the ordinary harmonic oscillator of frequency 
uj with a strong spin-orbit coupling term. This justifies the name of the Duffin- 
Kemmer-Petiau oscillator for our model as well. 



3 Exact solutions to the derived equation 

Because within the framework of the DKP equation the total angular momentum 
J = L + hs is conserved, the spatial variables for the components of the wave 
function H14f) are separated and we can write |22 

^=^^{r)YjMm, F=^-^fdr)Yf[,{n), 
G^lj29L{r)YfLm, K ^\Y.^dr)YfL{^), (17) 

where Yj[^{{l) are the so-called vector spherical harmonics. 

Substituting (|17|l into H15|l leads to ten radial equations which reduce to the two 
uncoupled sets associated with the (—1)'^ and (—1)''+^ parities. The (—1)'' solutions 
correspond to the natural-parity states. They are described by the functions (|17|) 
with the nonzero components fj,gj, kj±i. The orbital angular momentum for these 
states has a definite value J = L. The (—1)'^+^ solutions, with values of J = L — 1 
and J = L + 1 being mixed, are referred to as the unnatural-parity states. Their 
nonzero components are /,/±i, gj±i, kj, (j). 
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In the case of the natural-parity states, the radial equation has the oscillator-like 
form 

?i -nr H 5 y (bi - b2)hujTn 

-(61 + b,rmWr^ - ^^:lil±2l)fj^r) = 0, (18) 



whereas for the unnatural-parity states one gets a system of two coupled equations 

,d^ 
dr^ 



[n, + 5 6(bi + b2jnujm - [bi + 02) m uj r 5 )0(r) 



{n^i-^ + P— - {bi ~ b2)hcum - (bi - 62)"m"c^"r" - " '\ ' ^' )kj{r) 

ar c r 



E'-m^c^ ,^ ,,2222 ?iV(J+l), 



2bi^JiJ + l)t^E 
- 2 0r =0. (19 

Now we consider exact solutions to the last equations. It is easily seen that the 
obtained system decouples into two oscillator-like equations and, hence, has exact 
solutions if and only if 61 = or 62 = 0. For these two cases we have: 

(i) 61 = 1, b2 — (the value bi = const 7^ 1 can be absorbed into uj providing 
bi = 1). This choice of coefficients reproduces the DKP oscillator model proposed 
in Ref. 16. According to 118|l and l|19|l . its eigenvalues take the form 

E^-m^c* ^ {An + 2J + 2)nujmc^{L = J), 



E'^ - jr?c^ -f huj^m?d^ + 4J(J + \)E'^ = (4n + 2 J + 3);iwmc2(L J ± 1), 
E"^ -m^c" ~hiJ:^f^^^^^fT^J{TTV)E^ = (4ra + 2 J + 3);iwmc2(L = J ± lO?0) 

where n is the radial quantum number. Here the first equation describes the energies 
of the natural-parity states and the rest equations correspond to the unnatural- 
parity states. 

(ii) &i = 0, 62 = 1. With these coefficients we arrive at the alternative model 
for the DKP oscillator. Its eigenvalues are given by 

E'^-m^c^ = {4:n + 2J + 4)hu;mc^{L = J), 

E^-m^c^ = (4n + 2J + 6)?iwmc2(L = J±l), 

E^-m^c^ = {4:n + 2J + 2)hu;mc^{L ^ J ±1) (21) 

that can be rewritten with the common formula 

E^ -m^c^ = {4:n + 2L + 'i)hLumc^ (L = J - 1, J, ./ + 1). (22) 

In contrast to the equations (|20|l . the eigenspectrum of this model proved to be 
completely equidistant. Besides, these models differ in the zero-point energy. 

Note, that in the case J = 0, the radial equations should be deduced individually 
with putting f-i{r) = g-i{r) = fc-i(r) = into p7|l . For all this the expression 
for the eigenvalues with L = J + 1 = 1 remains still valid. 



4 Discussion 

Thus, we have considered the DKP equation with the interaction introduced through 
the non-minimal substitution of general form. This allows us to find a new exactly 
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solvable model for the spin-1 oscillator that possesses the completely equidistant 
eigenspectrum. 

Moreover, the eigenvalues of the proposed model differ from those of the spin-0 
oscillator [111 [T^ in appearance only the constant addition term Ahui. It means 
that within the framework of the alternative model the spin-1 particle behaves like 
the spinless one. In fact, as soon as we insert 61 = into the terms with the 
spin-orbit coupling vanish there. 

Hence, it appears that the exactly solvable models for the DKP oscillator con- 
sidered in the previous section differ dramatically from each other. This should be 
regarded as arising from the different Lorentz structure of the oscillator coupling 
terms in these models. 

For elucidating this difference let us turn to the Lorentz covariant form of the 
DKP equation ifT^ 

[c/3>M + ^6i(5^./?A + PxS^,,)T^''^ - lb2f3^,V5V'' - mc^]t^ = 0, (23) 

where 6*^1, = i{l3^l3i, — (3^13^) is the tensor of spin, and external potentials are 
involved in the Lorentz-tensor, T'^'^^, and Lorentz- vector, 1/^, interactions in spinor 
space. Then the equation H13|) is immediately restored with setting 

yO _ 2^000 _ jiijo _ rpiiui _ Q "\/* — — _rpiao _ ^^^^^i ^24) 

Taking into account (|23ll . we conclude that the model of the DKP oscillator 
with 61 = 1, 62 = considered in the paper |16| realizes the case of the Lorentz- 
tensor coupling, whereas the alternative model with 61 = 0, 62 = 1 proposed here 
includes the coupling of the Lorentz-vector type. 

A few remarks concerning the lack of the spin-orbit coupling are in order. There 
is another way to introduce the harmonic oscillator in relativistic quantum mechan- 
ics. It has been realized within the framework of the Dirac equation by means 
of mixing vector and scalar harmonic potentials with equal magnitude and equal 
or opposite sign (for references see Ref. 10), that results in a quadratic equation 
for each spinor component. The energy eigenvalues obtained in the such approach 
are characterized by the (2?! -f- 1) degeneracy and the lack of contribution from the 
spin-orbit interaction just as in our model. In addition, the lack of the spin-orbit 
coupling occurs also in the description of a spin-1 particle with the gyromagnetic 
ratio 5=1 moving in the external electromagnetic field . 
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